Nonanalytic corrections to the specific heat of a three-dimensional Fermi liquid 
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We revisit the issue of the leading nonanalytic corrections to the temperature dependence of the 
specific heat coefficient, 7(T) = C(T)/T, for a system of interacting fermions in three dimensions. 
We show that the leading temperature dependence of the specific heat coefficient j(T) — 7(0) oc 
T 3 In T comes from two physically distinct processes. The first process involves a thermal excitation 
of a single particle-hole pair, whose components interact via a nonanalytic dynamic vertex. The 
second process involves an excitation of three particle-hole pairs which interact via the analytic static 
fixed-point vertex. We show that the single-pair contribution is expressed via the backscattering 
amplitude of quasiparticles at the Fermi surface. The three-pair contribution does not have a 
simple expression in terms of scattering in particular directions. We clarify the relation between 
these results and previous literature on both 3D and 2D systems, and discuss the relation between 
the nonanalyticities in 7 and those in spin susceptibilities. 

PACS numbers: 71.10.Ay,71.10.Pm 



INTRODUCTION 



The thermodynamic properties of itinerant fermionic 
systems is a subject of long-standing experimental ^1 
3l and theoretical i[ 
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interest 
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It is generally accepted |3lL 




low-temperature behavior of a wide class of interacting 
fermionic systems is controlled by the Fermi-liquid fixed 
point. This implies that the leading temperature de- 
pendences are the same as for free fermions, but with 
renormalized parameters. However, the first sublead- 
ing corrections may differ qualitatively from those of 
non-interacting fermions. This phenomenon was first 
noticed in the context of the specific heat coefficient, 
7(T) = C{T)/T. For non-interacting fermions, 7(T) 
has a regular expansion in powers of T 2 about T = 0, 
so the leading temperature dependence is 7(T) — j(T — 
0) ~ T 2 . For interacting fermions, j(T) is not an an- 
alytic function of T 2 ; the leading temperature depen- 
dence is instead proportional to T 2 ln(T) in three di- 
mensions (3D) and T in 2D The T 2 lnT term was first 
found by Eliashberg in a theoretical study of electrons 
interacting with acoustic phonons Q , and the possibility 
of a nonanalytic temperature dependence of 7 was sub- 
sequently (but apparently independently) inferred from 
measurements of 7 for 3 He by Abel, Wheately and An- 
dersen [ij. The 3 He measurements led to a l arge litera- 

Efl l2ima ^ ^^jjiiM il? 1^[ilf il][which we review 
in detail below. We note here a few crucial (and seem- 
ingly contradictory) results, (i) In a very important pa- 
per [15j , Pethick and Carneiro showed that for a 3D Fermi 
liquid the T 2 InT term in 7 was a Fermi-liquid effect, as- 
sociated with a combination of the multiple scattering 



of particle-hole pairs with small total momentum and a 
particular behavior of the quasiparticle interaction func- 
tion / PiP +q (also at small q). (ii) More recently, it has 
been demonstrated [13] that in 2D, the prefactor of the T 
term in j(T)—j(T — 0) ~ T is determined entirely by the 
squares of charge and spin components of the backscatter- 
ing amplitude. [By "backscattering" we mean scattering 
of fermions with almost opposite momenta. Transferred 
momentum can be either small or near 2kp.] On the other 
hand the previous 3D work [lil llfil Hi| found no special 
role of backscattering for the T 2 In T nonanalyticity. ( Hi ) 
The older work left the impression that in 3D the nonan- 
alyticities were particular to 7 and did not contribute to 
susceptibilities, whereas more recent work (based mainly 
on perturbative calculations) has demonstrated that non- 
analytic corrections to the spin susceptibilit y o ccur both 
in 3D H3 and in 2D [13, HH US E3, [2^, [2^, [3^| . 

In this paper we clarify the relation between nonan- 
alyticities in 7 for 3D and 2D systems and make a few 
remarks concerning the relation between the nonanalyt- 
icities in 7 and those in susceptibilities. We demonstrate 
that backscattering plays a special role also in 3D, in the 
sense that a part of the T 2 In T term comes entirely from 
backscattering. This backscattering contribution evolves 
smoothly between 2D and 3D, and is entirely responsi- 
ble for the nonanalytic part of 7 for D < 3. In D > 3, 
however, there exists another, physically distinct, con- 
tribution to the T 2 lnT term in 7. This contribution 
does not occur for D < 3 and is not expressible solely 
in terms of the backscattering amplitude. We argue that 
this contribution is important for the 3D spin suscepti- 
bility as well. We also clarify the relation between the 
expressions for the prefactor of the T 2 In T term via the 
forward-scattering interaction in the Pcthick-Carneiro 
approach (and in subsequent analyses based on bosoniza- 



tion [22l |23|) and via the backscattering amplitude. In 
particular, Pethick and Carnerio expressed the nonana- 
lyticities in terms of the small q form of the scattering am- 
plitude /p iP + q for momenta slightly displaced from the 
Fermi surface. As we will demonstrate explicitly, their 
result can be re-expressed in terms of a dynamical in- 
teraction between quasiparticles at the Fermi surface, in 
which backscattering plays a crucial role. At the same 
time, we show, in disagreement with Refs. 0, 0] that 
processes involving forward-scattering between quasipar- 
ticles at the Fermi surface do not contribute to the non- 
analyticity in 7 although forward scattering does make a 
nonanalytic contribution to the self energy. 

The remainder of this paper is organized as follows. 
In section II, we outline the qualitative physics under- 
lying the results presented here. In section III, we 
present the formalism and introduce forward-scattering 
and backscattering amplitudes. Sections IV and V give 
detailed calculations of the nonanalyticity in the entropy 
for the cases of three and two spatial dimensions, respec- 
tively. Section VI presents presents a brief discussion of 
the spin susceptibility. Section VII explicates the relation 
of our results to previous calculations of the self-energy, 
and section VIII compares our results for the entropy to 
those obtained in prior work. Section IX presents a brief 
comparison of our results to the specific heat data for 
3 He. Section X presents summary and conclusions. An 
explicit calculation of the 2&f contribution to the non- 
analyticity in the specific heat is presented in Appendix 
A. 



QUALITATIVE PHYSICS 

In this section we present a qualitative discussion of 
the physics underpinning our key results. Subsequent 
sections provide the formal justification. To compute the 
low-T specific heat of a Fermi Liquid, one expands the 
thermodynamic potential (EE) in the number of particle- 
hole pairs excited above the Fermi-liquid ground state. 
Each excited quasiparticle interacts with the correspond- 
ing excited quasihole and with other particle-hole pairs 
via interaction vertices, T. Some low order diagrams are 
shown schematically in Fig. ^ 

To obtain the entropy S = —cE/dT we must average 
over all excited particle-hole pairs, weighting the contri- 
butions by the derivative of the Bose distribution func- 
tion driB (e) / dT, which is odd in e and, at low tem- 
perature T ensures that all excited electrons are near 
the Fermi surface. The contributions also involve the 
spectral function of the quasiparticlc-quasihole pairs [the 
imaginary part of the particle- hole propagator U(s,q)}, 
which is also odd in e, and, critically, is not an ana- 
lytic function of e. For example, in D = 3, ImII(e, q) ~ 
(e/vpq)Q(vpq/e—l) with 9 (x) the step function. There- 
fore, only states with an odd number of excited pairs 
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FIG. 1: Skeleton diagrams for the thermodynamic potential. 
The solid lines represent electron propagators, and the dashed 
lines represent the interaction potential. 

contribute to the entropy, and thus to j(T) — dS/dT. 
The leading term has only one excited particle-hole pair 
(diagram la in Fig^l and for momentum-independent 
interaction T yields for the interaction-dependent part of 
the entropy 

^^Tjdes d -^lJ d X (1) 

The momentum integral in Eq [Q is dominated by large 
q ~ kp, and so yields a constant, whereas the frequency 
integral yields 5S(T) cx T, i.e., 7 =const. This is the 
leading Fermi- liquid result. 

The subleading terms contain contributions from 
three, five, etc thermally excited particle-hole pairs. For 
thermally allowed excitations e ~ T , so each extra 
particle-hole pair brings an extra factor of T to the ther- 
modynamic potential. One might therefore expect an ex- 
pansion of the form 8S = jqT + aT 3 +... . However, the 
non-analytic behavior of the particle-hole spectral func- 
tion implies that this conclusion is not correct. The non- 
analyticity enters in two ways. First, processes involving 
three pairs with the same total momentum q (see, e.g., 
Diagram 3c of Fig^l contain the cube of the spectral 
function, so give a term proportional to (e/q) 3 ; the in- 
tegral over q leads to a logarithmic divergence cut off at 
e ~ T , and hence, in 3D, to a T 3 In T contribution SS(T). 
Second, the general form of the interaction vertex T has 
a term of order (e/q) 2 arising from the Kramers-Kronig 
transform of ImII(e, q) which, in combination with the 
e/q from a single excited particle- hole pair, produces in 
3D an additional T 3 InT contribution to SS(T) (see, e.g., 
diagram 2b in Fig. QJ. 

Note that dimension D — 3 is marginal in the sense 
that only for D < 3 the leading nonanalytic term is larger 
than the leading analytic <D(T 3 ) correction. As will be 
seen, it is also marginal in the sense that for D < 3 only 
the single-pair term contributes to the nonanalyticity. 
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Both contributions were identified by Pethick and 
Carneiro . They can be also clearly seen in the "para- 
magnon model" [3 of electrons coupled to over- 
damped spin fluctuations (although only the three-pair 
contributions were discussed in [lOl UJ). The param- 
agnon part of the entropy is obtained by summing up 
loop dia gram s with a small momentum transfer and is 
given by 



gran 



6S = 



ds 



8irT 2 J sinh 2 (e/2T) 



(27T) 



D 



tan 



2ghnIl(E,q) 
1 + 2 5 Rcn(e, q) ' 



(2) 



where g is the coupling constant in the paramagnon 
model (see Sec. D below). One obtains one (e/q) 3 term 
from the process involving three real excited particle hole 
pairs by treating ReLT as a constant and expanding the 
tan -1 to third order in Imll ~ e/q. In addition, one 
obtains another (e/q) term by combining one power of 
Imll cx e/q (representing a real excited particle-hole pair) 
with the nonanalytic term in Rell oc (e/q) 2 (representing 
the nontrivial dynamic structure arising from ReLT). 

The two processes evidently represent different physics 
and differ mathematically as well. The single-pair mech- 
anism occurs already at the second-order in the interac- 
tion g and is not specific to D = 3. A simple extension 
of the 3D analysis to arbitrary integer D shows that the 
single-pair mechanism yields a non-analytic 5S cx T D in 
any dimension (with an additional InT in odd dimen- 
sions). On the contrary, the three-pair mechanism oc- 
curs to third order in the interaction and gives rise to 
non-analyticity in SS only in odd dimensions D > 3. 

Let us consider the kinematics of these contributions in 
more detail. In both processes, the nonanalyticity arises 
from Landau damping, which for the physically relevant 
limit of e <C q is dominated by fermions with momenta 
nearly perpendicular to q. For a single-pair process, rep- 
resented, e.g., by diagram (Fig.n2b), the constraint that 
the momenta on all four fermion lines are perpendicular 
to the same vector q implies that in two dimensions all 
four momenta are either nearly parallel to each other or 
two are nearly antiparallel to the other two. A closer 
analysis [13] shows that only the nearly antiparallel case 
contributes to the nonanalyticity in 7 so the nonanalyt- 
icity is controlled by backscattering. However, it appears 
surprising that backscattering would play a special role 
in 3D. Indeed, in D = 3 the constraint that all four mo- 
menta are perpendicular to q implies only that all four 
momenta lie in a common plane, so that any multiparticle 
excitation would seem to involve a more general scatter- 
ing process involving four fermions, two of which have 
momenta near k and two near p but with k and p per- 
pendicular to q but not otherwise constrained. Indeed 
the three pair contributions have precisely this structure. 
Nevertheless we will show explicitly that in the single pair 



process only the case of nearly antiparallel momenta con- 
tributes. 

Some physical understanding of the importance of 
backscattering even in 31? may be obtained by consid- 
ering the problem in a slightly more general way. The 
discussion just given has focussed on "small momenta" 
in the sense that it considered processes describable in 
terms of electron-hole pairs where the electron and the 
hole momenta are very close (on the scale set by kp), 
and the nonanalyticity arises from the singular structure 
of the Landau damping ~ e/q in the small e, q limit. 
However, Fermi liquids also exhibit nonanalyticities in- 
volving momenta near 2kp , (reflected for example in the 
slow decay of Friedel oscillations), and previous literature 
has raised the question of the contribution of 2k p pro- 
cesses to the nonanalytic behavior [l^. To understand 
this issue it is useful to consider again the simple per- 
turbative contribution shown in diagram 26 of Fig. ^ 
which may be written as the product of two particle-hole 
bubbles: 



26 



ds 



dn B (e) 
dT 



d D qlmU(q,e)ReU(q,e). (3) 



There are two possible momentum regions which may 
give rise to nonanalyticities: small momentum and q ~ 
2kp. Consider the contribution to f^b arising from large 
momentum transfer, q = (2k f + q)q. The nonanalyt- 
icity arises from the "Landau-damping" structure of n, 
which in 3D and for momenta near 2kp is Imn(q, e) cx 
e 6(|e| - q) and Ren - Ren(2/c F ,0) - e 2 /q (see Ap- 
pendix A). We note that the orientation of q is irrele- 
vant, so J d D q becomes a one-dimensional integral over 
the scalar quantity q leaving a logarithm for D = 3.. 
The 2kp "Landau damping" terms involve fermions with 
momenta close to kp, so that if q ~ 2k f then in each 
bubble one has one fermion of momentum k w q/2 and 
one with k « —q/2; thus the 2kF contribution involves a 
one-dimensional process controlled by the "backscatter- 
ing amplitude" T(k, — k; k, — k). However, we may also 
view the diagram for f^fc in a different way, regarding the 
two fermion lines with momentum k ~ q/2 as a bubble 
with a small total momentum q, and similarly with the 
two lines of momentum near — q. Thus for the process 
with one real excitation and first dynamical correction to 
interaction, the 2kp contribution may be subsumed into 
the small q contribution. 

The total nonanalytic contribution to SS(T) may thus 
be described in terms of two sorts of small momentum 
processes. One may be thought of in terms of 2kF pro- 
cesses, and involves the vertex T(k, — k; — k, k). The 
other arises from consideration of small q processes and 
seemingly involves the effective interaction T(k, p;k, p) 
with arbitrary values of k • p. Graphically, the two con- 
tributions differ only by interchanging the two outgoing 
momenta in the effective interaction. It is thus plausible 
to expect that both contributions are expressed in terms 
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of the spin and charge components of the same antisym- 
metrized vertex. The vertex r(k, — k; — k, k) is the spin 
part of the backscattering amplitude. We may expect 
that r(k, p; k, p) is expressed in terms of charge compo- 
nent of the same amplitude. If this is the case, then p 
must be antiparallel to k, i.e., the original small q, single- 
pair contribution also involves ID process, and hence the 
total single-pair contribution involves ID scattering. 

For the three-pair mechanism, this argument does not 
hold. The kinematics of 2kp processes is such that if 
more than two pairs are involved the result does not have 
a strong enough nonanalyticity. Indeed, combining three 
factors of ImiT(£, 2kp) cx e, we obtain only a regular term 
in SS(T). Thus the argument for the presence, in the cal- 
culation, of a purely one-dimensional scattering process 
fails and backscattering plays no special role. 

A principal aim of the present paper is to present a 
more rigorous analysis substantiating the qualitative ar- 
guments given above. In this analysis the effects of single- 
pair and three-pair processes are separated and the con- 
tribution of each to the nonanalyticities in 7 is deter- 
mined in two and three dimensions. 

FORMALISM 

We calculate the specific heat coefficient j(T) = C/T 
by evaluating the difference between the thermodynamic 
potential, 3, and its value at T = 0. The Luttinger-Ward 
expression for 3 is |33| 

S = -Tr(ln[G - 1 -S] + EG) + 3 skol (4) 

with E the exact self-energy (regarded as a functional of 
the full Green function G) and S s k e i is the usual skeleton 
diagram expansion for all interaction corrections to 5 
that are not accounted for by the first two terms. The 
trace is taken over space, time and spin variables. 

It is not necessary to obtain the nonanalytic contri- 
butions the self-energy E in order to evaluate the non- 
analytic contributions to 3. For what appear to be his- 
torical reasons associated with the fact that nonanalytic 
contributions to £ were a focus of many earlier studies, 
the previous literature evaluates the Trln [Gq 1 — E] sep- 
arately from the other contributions, thereby effectively 
adding and subtracting contributions arising from non- 
analyticities in £. However, the stationary of 3 with 
respect to variations of E implies that it is in fact suffi- 
cient to consider only the skeleton diagram 5 s k e i, which 
may be evaluated using the leading renormalized quasi- 
particle Green function G qp = Z^ 1 [u — v* F (k — kp)] 1 , 
where Z is the renormalization factor and v* F is the renor- 
malized Fermi velocity. We emphasize that, in physical 
terms, this corresponds to expanding 3 in the number 
of particle-hole pairs thermally excited above the ground 
state, but assuming that these interact according to the 
T = Landau Fermi Liquid fixed-point Hamiltonian. 



A further subtlety arises: even if the Green function 
takes the fixed-point form, the 2kp particle-hole bubble 
will have a nonanalytic temperature dependence. How- 
ever, as we discussed in the previous section, singular 
2k p scattering can be re-expressed as small momentum 
scattering, and is fully accounted for if we restrict to 
small q but consider wavy lines in Fig. ^ as antisym- 
metrized interactions, i.e., as fixed-point vertices. The 
small-g diagrams have nonanalyticities arising from the 
Landau damping e/vpq which has negligible temperature 
dependence. Once the restriction to small q is accom- 
plished, we may calculate the leading low-T behavior of 
3(T) — 3(0) by focusing only on the T-dependence aris- 
ing from the difference between T and J cko /2ir, and 
neglecting any explicit temperature dependence of G, E 
or interaction vertices. These latter contributions give 
only analytic, O (T 2 ) corrections to 7. Appendix A con- 
firms this argument by presenting an explicit calculation 
of the contribution to S arising from 2kp processes. 

One can employ two strategies to calculate nonanalytic 
terms in 3 (T) . The first one is to evaluate the ther- 
modynamic potential directly in Matsubara frequencies 
e m . This strategy was adopted in recent study of 2D 
systemsj^. A drawback of this approach is that it docs 
not distinguish between the real and imaginary parts of 
U(e,q), and the physical picture of excited particle- hole 
pairs does not arise. The second -and more intuitive- 
approach, adopted in this paper, is to work in real fre- 
quencies, when Imn oc e and Ren cx e 2 have different 
physical meaning: the former describes real particle-hole 
pairs with given momentum and energy, whereas the lat- 
ter describes the interaction arising from virtual pairs. 

To evaluate non-analytic terms in 3 for a generic Fermi 
Liquid we need two quantities. The first one is the 
propagator "P p h(£, q', n^), describing a particle-hole pair 
with small total momentum q and small energy e, re- 
spectively, and with the particle and hole momenta near 
the Fermi surface and in direction rife. The second one 
is the fully renormalized, particle-hole reducible vertex 
r Q , / 3 ;7 a(e, g; rifc, n p ) describing scattering of one particle- 
hole pair state into another. 

For small q and s, V p h may be written as a sum of two 
terms: "P p h = Vq + V. The analytic part Vq is a function 
of e and q 2 and contains contributions from virtual pro- 
cesses involving states both near and far from the Fermi 
level. The nonanalytic part, denoted by V, is determined 
solely by the properties of Fermi surface states and de- 
pends only on z = v* F q/e. Expanding the product of 
two quasiparticle propagators near the Fermi surface, we 
obtain 

p(z;n k ) = -±- Z (5) 

Sd z 1 - n fc • n q 

with Sd being the surface area of the sphere of unit radius 
in D dimensions. We absorb the factor of inverse velocity 
and the factors of quasiparticle weight into the definition 
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of r, and we have specialized to rotationally invariant 
systems. The frequency-dependent part of the polariza- 
tion bubble n(e, q) is obtained by averaging r P(e,q;n k ) 
over n^: 



U(z) = / dn k V{z; n fc ). 



In 3D, 



z + 1 



Ren(z) = —In 



= --e(|*|-l). 



(6) 
(7) 



In the perturbation theory, when the interaction de- 
pends only on the momentum transfer but not on the 
incoming momenta, one needs to know only the angle- 
averaged bubble. However, the interaction in a generic 
Fermi Liquid is described by a vertex T a p. n $(£, q; n^, n p ), 
which depends not only on the momentum transfer q but 
also on the relative directions of the incoming momenta 
ilk and n p . Therefore we need need the n^-dependent 
propagator V(e, q; nfc). In evaluating 5, we will then have 
to perform angular integrations of the convolutions of V 
and T. 

The fixed-point vertex r Q ^ ;7< 5(e, q\ n^., n p ) is an 
(anti)symmetrized sum of all scattering processes be- 
tween fermions at the Fermi surface. We may neglect 
any analytic dependence of T on s and q separately, and 
consider T to be a function only of z. The fixed point 
vertex is also a tensor in the spin space. It is conve- 
nient to express this tensor as a sum of charge and spin 
components 

T a p- n s(z\ n fc , np) = 5 ai 8ps^c{z] n fc , n p ) 

+Cq7 • CT /35 r s (z;n fc ,n p ), (8) 

where c and s refer to the charge and spin sectors, respec- 
tively. In a model with bare fermion-fermion interaction 
U{q)p{q)p{—q) to first order in the interaction we have 

r Q /3 ;7 5(^; n fc , Up) = S aJ 6i35u(Q) - SusS^U^kFlllk - n p |) 

(9) 

and thus 

F c (z;n fc ,n p ) - u(0) - (l/2)u(fe F |n fe - n p |) (10) 
r s (z;n fc ,n p ) = -(l/2)u(fc F |n fc - n p |), (11) 

where u(q) = U(q)k F /(ir*v%). 

Pethick and Carneiro [13 showed that the charge and 
spin components of F contribute independently to the 
nonanalyticity in S, i.e., Sna = ^2 a w a S , where a = c, s 
and w c = 1, w s = 3 (for SU (2) fermions). We explicitly 
verified and confirmed their result, which indeed follows 
simply from an elementary consideration of the possible 
distributions of charge and spin vertices along fermionic 
loops. To make our presentation more compact, we will 





FIG. 2: Graphical representation of Eq. ill'Jt . 



consider only the charge component of F and omit the 
index c. The spin component will be restored in the final 
results. 

The forward-scattering and backscattering processes 
in these notations correspond to r(z;n p ,n p ) and 
T(z, n p , — n p ), respectively. In the first process, all four 
fermionic momenta are almost equal; in the second one, 
the two incoming and the two outgoing momenta are al- 
most antiparallel. 

The vertex r(z;n p ,rife) satisfies the integral equation 

r(z;n fe ,n p ) = r fc (n fc • n p ) 

+ J dmT k (ii k •n i )P(2;n i )r(z;n / ,n p ), (12) 

where the integral is over the area of the unit sphere 
in D dimensions. Virtual processes which contribute to 
Vo are absorbed into the "bare" vertex T k (iik ■ n p ) = 
r(oo;n/,n p ). The notation T k , borrowed from Ref.Q, 
comes from the observation that V vanishes at z = 00, 
i.e., at s/vpq — ► 0. We recall that T k coincides, up to the 
Z— factor, with the quasiparticle scattering amplitude. 
The scattering amplitude differs from the Landau func- 
tion, which is related to F" (n^ • n p ) = T(z — 0, n^, n p ). 
The scattering amplitude may be expanded in spherical 
(or, in two dimensions, circular) harmonics as 

r fc (n fc ,n p ) = T k {n k ■ n p ) = ^F L P L (n fc • n p ). (13) 

L 

Here Pl are Legendre polynomials in D = 3 and cosines 
in D = 2. The partial amplitudes, Fl, of T fe are related 
to the partial amplitudes, T^, of via 

r. 



l+r £ (2L + l)-i 



and 



i + r (2-^, )- 1 

in 3D and 2D, respectively. 

To generate an expansion in a number of particle-hole 
pairs we need an auxiliary vertex, T(z; n kl n p ), which sat- 
isfies the same equation as in Eq. l|12|l but with Re'P(z; n/) 
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Here the factors of Imp represent the excited particle- 
hole pairs and the factors of A (determined by combi- 
nations of the f from Eq. I12() describe the interaction 
between them. We note that because ImP is odd in e 
whereas T is even, only terms involving odd powers of 
ImP, i.e., odd number of pairs, contribute to SSna- To 
extract the nonanalytic behavior it is convenient to in- 
troduce x = e/2T and rewrite Eq. Ijl5|l in dimension D 
as 

2 D+1 S D T D f°° dx 



5Sna 



Vpk F /Tx 



„D+1 



7r sinh x 



(17) 



Here, we used the fact that lvaP{z, n p ) = if \z\ < 1. 



FIG. 3: Diagram 2b for the thermodynamic potential with 
fully dressed vertices T(q, \bfk — p), instead of the bare inter- 
action U(q) 



instead of full P(z; n;) 

r(z:n k ,n p ) = T k (n k ■ n p ) 

+ J dmT k (n k ■n / )ReP(z;n ; )f(z;n / ,n p ). (14) 

The leading low-T contributions to the entropy in- 
volve particle-hole pairs excited above the Fermi sur- 
face and can be mathematically described by combina- 
tions of ImV, representing the excited particle- hole pairs, 
and vertices T, describing the interactions between these 
pairs. The diagrammatic expansion of the thermody- 
namic potential in series of particle-hole pairs is obtained 
from the skeleton diagrams of Fig. ^ by replacing the 
wavy lines by fully dressed interaction vertices. For ex- 
ample, diagram 2b of Fig. \T\is replaced by the diagram 
in Fig. |21 An important new feature of this expansion 
is that the vertices now depend not only on the mo- 
mentum transfer (q) but also on the incoming momenta 
(k and p). Consideration of diagrams involving exci- 
tations and multiple scattering of pairs with total en- 
ergy e and momentum q leads to an expression for the 
nonanalytic contribution to the entropy per unit volume 
8S NA = -(l/V)dS ske i/dT of the form 



where 



de f d D q 



7T J (2tt) d 4T 2 sinh 2 

C_iy+i 

$ (£jg ) = <i>( z) = ]rL_z_ x 

f[ J dn p ImP(z,n Pj )A(z,{n P:j }) 



$(£,<?), (15) 



(16) 



NON-ANALYTICITY IN SPECIFIC HEAT IN 
DIMENSION D=3 

Overview 

This section presents results for the nonanalytic en- 
tropy contribution <5Sna oc T 3 In T for a three- dimen- 
sional Fermi Liquid. After an overview we present re- 
sults based on expanding Eas. 1)121 117|l in the number of 
physical excited particle-hole pairs. We show that (SS'na 
can be expressed in terms of series of partial amplitudes 
T^; although a compact, useful closed- form expression 
exists only for a part of <5Sna corresponding to a single- 
pair mechanism. In the last subsection we consider a toy 
model in which the expansion in particle-hole pairs can 
be carried out to all orders. 

The key point of the calculation is this: the require- 
ment of odd powers means that in Ea. ll6|) only terms 
with odd j, i.e., odd powers of ImV(z,n Pj ) contribute 
to &(z) . At large z, ImP is of order z~ l ; while 
A = A + Bz- 2 + 0(z- 4 ) is finite at z" 1 = 0. The leading 
large- z contribution is thus O the integral in Eq. 

(|17fl is dominated by the upper limit and gives a contri- 
bution to the 2-integral ~ T~ 2 +const; the net result is 
contributions of order T and T 3 to the entropy which are 
not of interest here. The subleading contribution come 
from the j = 3 term in Ea. H16J) and from the z~ 2 term 
in A combined with the j = 1 term in Eq \Fg\ . Both give 
contributions of order z~ 3 to $(z). 

Expanding <fr(z) asymptotically for large argument 
yields $ (z) at large z as 



<f>(z -> oo) = — + ^ + (18) 

z z A 

Subtracting off the leading term, inserting the result in 
Eq. (|T7ll . and integrating gives for the entropy per unit 
volume 

3 



ss 



4tt 3 



N A 



T 



In 



T 



(19) 
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The problem is therefore to calculate $3 . It is convenient 
to split $3 into contributions from processes with one and 
with three excited real particle- hole pairs, i.e., 



$ 3 = $f ) 4 <I» 



(3). 
3 ! 



and to discuss the contributions separately. 



Expansion in the number of excited particle-hole 
pairs 

Single particle-hole pair 

For the term involving only one excited particle-hole 
pair, reference to Fig. [T]p, shows that the interaction term 
is the quasiparticle vertex f , from which we need the 
contribution proportional to z~ 2 . Hence 



z 3 f 

$y = iirn^— / dn p lm.V(z;n p )r(z;np,n p ) 



(20) 



with T(z: n p ) being the order z~ 2 term from the solution 
of Eq. H12JI . Direct inspection shows that there are two 
contributions to T(z; n p ); one, denoted of first order 
in ReP and another one, denoted T^- 2 \ of second order. 
These contribution are give by 



dn k [T k (n p -n k )] 2 ReV(z,n k ), (21) 



T^(z) = J dn kl dn k2 T k (n p ■ n kl )ReV(z,n kl ) 
xT k (n kl -n k2 )ReV(z,n k2 )T k (n k2 -n p ). (22) 

Accordingly, $3 is a sum of two contributions, $3 
and $ 3 1,2) . If T fe is angle- independent, the angular inte- 
grations in Eqs.|2TE2l yield T^(z) ~ (r fc ) 2 ReII(z) cx 
z- 2 and T< 2 )(z) - (r fc ) 3 [Ren(z)] 2 cx z" 4 , respectively 



[cf. Eq.®]. Hence $ 



(1,1) 



product n/j - rip explicitly in terms of polar angles referring 
to the direction of g, i.e., 

T k (n k -n p ) = 

f „ (cos 9 k cos 6 p + cos <fi k sin k sin 9 p ) n . (23) 

n 

In these notations, the backscattering amplitude is 

r BS - r fe (n p , (-n p )) = £(-l) n f n . (24) 

n 

Inserting Eq. ill' Ml) into Eq. I|21|) and then the result into 
Eq. J2JJl, recalling that ImP = -z^ 1 5{z^ 1 -xtp -Ta. q )/A-K, 
exploiting the azimuthal symmetry and the reflection 
symmetry about the plane n p • n q — 0, and defining 
y = cos6» fc , yields 



(1,1) _ 



■- Vf f 1 



where 



f 1 z' 1 

I nm = lim z 2 / dyJ nm (z,y)P — : 

z^oo \ x z — y 



and 

Jnm \Z- y) 



f 7T (- + cos 0\/ 1 - z 2 \/l- y 2 

J 2n \z 



(25) 



(26) 



L—Q...n+m 

X 



(n + m - L)\L\ 



(zy [vi^^v^y^ 1 L . (27) 



0(1) and % ' ' tx z- 2 < 
However, if T k is angle-dependent, <&z~'^ and 
$3 ' are of the same order. Indeed, the large- z limit 
of Re"P(z, rife) in Eq.JSJ) is ~ z^ 1 / (n k ■ n q ), which means 
that the product of two ReP(z, rife) in Ea. (|22|l is odd in 
n kl 2 . However, this oddness is compensated by a com- 
bination of angular harmonics of T k in Ea. l|22JI . which 
contains a factor of (n^ • n 9 )(rifc 2 • n 9 ). As a result, the 
angular integrals are finite and r^ 2 )(z) cx (ReP) 2 cx z~ 2 , 
which is of the same order as rW(z). One can readily 
verify that terms with three, four, etc factors of ReP are 
irrelevant, as they scale at least as z~ 3 . 

We consider the two contributions, $3 and ^' 2 \ 
separately, beginning with the first one. A conventional 
way to proceed would be to expand r(rifc -n p ) in Legendre 
polynomials, as in (|13|) . and evaluate the integrals order 
by order in T^. However, it is more convenient to expand $( 1 > 1 ) 
r(rife • n p ) in powers of • n p , and to write out the scalar 



If n + m = 2P is even, then only terms with even L 
contribute to the integral in Eq.J57|), and J is an even 
function of y. In this case, the principal value integral 
gives a result of order z -2 , and only the L — term in Eq. 
(|27|l should be retained. If n + m — 2P + 1 is odd, then 
J is an odd function of y. In this case the principal value 
integral gives a result of order z _1 and only the L = 1 
term is needed. The (f> integral may then be performed, 
leaving 

_ (2P)\ f ,„ z{l-y 2 r 
l n +m=2P — n ,p, nh , / ayF— — j — (2») 



2 2P (P!) 2 
(2P+1)! 



z -y 

.2^P 



dy(l-yY. (29) 



1 n+m=(2P+l) 2 2P (P\) 2 

Finally, performing the y integrals leads to 

Inm = (-l)" +m - (30) 

Substituting this into Ea. (|2*5jl . we find 

2 



bs- (31) 
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We see that the contribution $3 is expressed solely 
in terms of the backscattering amplitude, i.e., it comes 
exclusively from ID scattering processes. Notice that 
Eq. (|31|) can be also expressed in terms of the angular 
harmonics of T k , introduced in Eq. (|13|) , as 



4> 



(1,1) 



1 BS — 



5 fe^y 



(32) 



We pause here to emphasize the non-triviality of this 
result. As we have already mentioned, the same process 
with one excited particle-hole pair leads to a nonanalyt- 
icity in D = 2. For two dimensions, however, an simple 
geometric argument, which we displayed in Section II, 
shows that ID processes determine the nonanalytic part 
of the entropy. The straightforward extension of this ar- 
gument to D > 2 suggests that all relative angles are 
important, but Ea. l|31[) shows that ^^'^ still comes only 
from ID scattering. 

We next consider the contribution involving I^ 2 ). Fol- 
lowing the arguments given above we find 



(1,2) 



4 ^2 

n.m.l 



Calculations along the same lines as above yield 
with 



'nm,!=(2P+l) 2 2P (P!) 2 " V2 



2P+1! (n \ /m \ 



I{m) = (l-y 2 )" 1 - 



(33) 



(34) 



(35) 



The integrals may be performed in terms of T-functions 
but the final expressions seem complicated and not very 

useful in general. For further comparison with Rcf. 151, 

(1 2) 

we present here the result for $3 , obtained by trun- 
cating the sum in Ea. (|33|l at n, m, I = 4 and writing the 
result in terms of the angular harmonics of T k : 



(1,2) 
3 

7T 

12 



7T 
4 



fir 2 , 



rr 2 2r 



^0X2 - 
1 8 - 

3flf 3 + yf 



-f 2 
• 1 2 



(36) 



Three particle-hole pairs 



The second contribution to the nonanalytic term in the 
specific heat of a three-dimensional Fermi Liquid comes 
from terms in Eq. I|19f) involving three excited particle- 
hole pairs ( e.g. diagram 3c in Fig. This contribution 
gives the "paramagnon" term in the specific heat stud- 
ied, e.g., in Ref. |lOj. In this case the quasiparticle in- 
teraction A involves three factors of the static scattering 



amplitude T k whereas ImP can be taken in the z — > 00 
limit, yielding 

$3 3) = i§2 / dn pidn P2 dn P3 S(n pi ■ n q )6(n p2 ■ n q ) 

xS(n P3 • rig) r fc (iV • n P2 )T k {n P2 • n P3 )r fe (iip3 . n pi ). 

Substituting Eq. ill'oli into the equation for $ 3 3 ^ and eval- 
uating the integrals, we obtain 



*<» = 



7T 3 



(I + m)\(l + n)\ 



Iran 



(37) 



Truncating the sum at l,m,n = 2 and re-expressing the 
result in terms of Tl [Eq. (|13Jl ] , we obtain 



(3) 



48 



Fo + i F 2 



1-, 27 -o 

-r? + — 

4 1 256 2 



(38) 



The final result for the nonanalytic part of the entropy in 3D 



Combining all three expressions $3 1 ' 1) ! $^ 2) and^ 3) , 
and restoring the sum of charge and spin components of 
r , we finally obtain for the entropy per particle 

3 



5S nA (T) = ~(^- 
5 \v F k F 



In 



v* F k F 
T 



(K c + 3K S ), 



K a — r 2 B S + r a,ir 2 ,0 + ^l,l^a,2 



F2 ( ( r ^° + I r ^ 2 



1-, 27 

_p3 J p3 

4 a ' x 256 °' 2 



.(39) 



where a = c, s, r a .BS = r a (oo; n^, — n^), and dots stand 
for the terms involving r Qj „ with n > 2. 

A comment is in order here. In the consideration above 
we assumed that the backscattering amplitudes r aj BS are 
temperature-independent. Strictly speaking, this is not 
the case as the backscattering amplitudes describe pro- 
cesses with total zero momentum and thus can be re- 
expressed via the partial components of the pairing ver- 
tex [34| • One of these partial components diverges at the 
pairing instability which, for a repulsive interaction, is 
of the Kohn-Luttinger type [35}), so our consideration is 
only valid at T > T c . Even at these temperatures, the 
backscattering amplitudes acquire a logarithmic temper- 
ature dependences from the Cooper channel, and behave 



as 



r Q: Bs(7 1 ) CX V] 7 , a 1 rr/F V 

rT^O ( a n + PnlnT/E F ) 



(40) 



where a n and (3 n are constants. At T c , a n +/3 n lriT/E F = 
for a particular n = uq. In the rest of the text, we ne- 
glect this complication and assume that the system is 
substantially far away from T c that both spin and charge 
components of the backscattering amplitudes can be ap- 
proximated by constants. 
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Ring Diagram (RPA) Approximation 

In this subsection we present results for a toy model 
in which the expansion in the number of excited particle- 
hole pairs can be carried to all orders, in order to demon- 
strate explicitly that only the squares and the cubes of 
the partial components of the scattering amplitude, but 
not higher powers, determine the non-analyticity in S(T). 

We consider an artificial model of fermions with spin 
degeneracy N — > oo, coupled by a contact interac- 
tion C/(n p - rife) = Uq + U\n p ■ rife. Standard tech- 
niques yield, for the large-N limit of the entropy Soo — 
-Hm N ^ 00 {dE{N)/dT)/N 



Sao — Sf 



-Im 



'ring 

de 



(41) 



7T AT 2 sinh ^ 



d 3 q 
(2^3 



x Trhi[(l-r(n p -n fc )7?(z,n fc ))], (42) 

where Sf — —d/dT [TrlnGp 1 ] the free fermion entropy, 
the trace taken over angle (Tr = rjdnfc), T = U{n p ■ 
n k )k 2 F /{n 2 v F ) =T Q + Yxn p ■ n k and 



V ph {z,n k ) = - 



1 cos(n fe -n <z ) 



(43) 



4-7T cos (rife ■ n g ) — z 1 ' 

Note that V p h is the full propagator; we do not split it 
into analytic and nonanalytic parts. Making the changes 
of variables to x = e/(2T) and z — v F q/e, and inte- 
grating over x as before we obtain for the entropy per 
particle 

2tt 3 / T \ 3 f v l Tx 



S, 



ring 



Z 



■dz 



5 \v F k F / ji 
x ImlnDet [1 -T(np ■ n k )V{z,n k )} . 

To evaluate the determinant we must find the eigen- 
values of the operator 1 — TV, i.e. we must solve 

AA(z, rip) = A{z, n p)~J dn k Y {rip ■ n fe ) V ph {z, n k )A{z, n fe ). 

(44) 

Writing 

A(n p ) = A + Ax cos8 p + A 2 sin0 p (e^ + e"*) (45) 

and using Ea. H23(l . we find that the determinant may 
be written as the product Det = D\ X D\ (ImlnDet = 
ImlnL>i + 2ImlnL» 2 ), with 

d 1 = (i + r /oo)(i + r 1 / 11 )-r r 1 / 2 1 (46) 

D 2 = l + IVn (47) 

and 

" x d(cos6») (cos0) n+m+1 



COS — Z," 



III 



1 d(cos9) cos6»(sin<9) 2 



(48) 
(49) 



Evaluating I mn , we find 
1 



Inn - 1 



, l-Z-' 7T 

— In- T- + i — 

2z 1 + z- 1 2z 



hi 



hi 



(50) 
(51) 
(52) 



3 z 2 w 

hi = {loo -hi)- (53) 

Substituting Eqs.J53Jl into Eas. (l46l47fl and simplifying 
yields 

Di = A (1 + f (Ioo - 1)) (1 + ri(/n -1/3))) 

-fofi/o 2 ! (54) 
D 2 = B + ri(Jn - 1/3) (55) 

1 + Ti/3, and the 



with a = (i + r )(i + r 1 /3), B 

reducible amplitudes Y given by 

r 



r 



fi 



i + r 

Yi 



1 



3 



(56) 
(57) 



The z integral in Eq |44j is dominated by its upper 
limit. The leading term is C(z _1 ) and gives a renormal- 
ization of the effective mass. The next term is (D{z~ 3 ) 
and gives the logarithm we require. Expanding in 1/z, 
we obtain 



Im In D\ 



Im In £>2 



inY iirYt in {Y 2 - 2Y Y 1 ) 



2z 3 



2z 

24z 3 ' _ 2z 3 
iirYi inYi 



2z 3 



mY 2 ^ 



nrY 2 



inY\ 



(58) 
(59) 



Az Az 3 Az 2 192z 3 ' 
Combining l|58l) and l|59l) to obtain ImlnDet = 
ImlnZ?i + 2Imln_D2, extracting the z~ 3 term and per- 
forming the integral in Eq. (|44(l yields for entropy per par- 
ticle 

3 



'ring 



TT 

~5~ 



In- 



vk F 



(60) 



To-rO + r 2 ri - 



12 



7T 2 f 3 

48 



This expression is in agreement with our previous result, 
Eq l|39|l once we note that to the order in N to which we 
work, only the charge component of Y contributes to a 
nonanalytic SS{T), and that in the RPA approximation, 
the dimensionless Y coincides with the Landau interac- 
tion function (this also implies that r„ are the partial 
components of the scattering amplitude). Note that, as 
before, the terms with coefficient it 2 in l|61|) arise from 
processes with three excited particle-hole pairs while the 
others involve processes with one excited particle-hole 
pair, along with a non-analyticity in the quasiparticle in- 
teraction. 
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NON-ANALYTICITY IN SPECIFIC HEAT IN 
DIMENSION D=2 

Overview 

For completeness, we here present results for the non- 
analytic contribution SS oc T 2 to the entropy a two- di- 
mensional Fermi Liquid. This contribution was found in 
Ref. |27J, |28j by exact evaluation of diagrams up to third 
order and by showing that results obtained by analy- 
sis of leading singular behavior of diagrams occurring in 
the general case matched the perturbative results exactly. 
The calculations were also performed in Matsubara fre- 
quencies. Here we show how the same result is obtained 
using an expansion in the number of real particle-hole 
pairs. In addition, we reconfirm the result that the non- 
analyticity in 2D involves only backscattering processes. 
Unlike the 3D case, this is the only non-analytic contri- 
bution to S(T). 

We begin from the two dimensional version of Eq 
[T7j in which the essential object is J^F kF / Tx zcIz<&2d{z). 
Again we expect the leading large z behavior of $ to 
be $2D ~ 1/z giving the Fermi- liquid behavior; but the 
different ultraviolet behavior of the momentum integral 
means that the remaining terms are ultraviolet conver- 
gent, so that the z integral is just a constant and the 
temperature prefactor gives the nonanalytic result in 2D: 
7 = d(SS NA )/dT ~ T. 

A further important question arises from the analytic 
structure of EL In D = 2, 



ReII(z) 
ImII(z) 



e(i-N) 

sgnz9 (\z\ 



1 



(61) 
(62) 



The singularity at the boundary of the particle-hole con- 
tinuum \z\ = 1 has be treated with care when expanding 
<5Sna in powers of ImEI, and indeed raises the possibil- 
ity that the z integral has appreciable contributions from 
the neighborhood of z = 1 , where the kinematics of Lan- 
dau damping does not guarantee that the contributions 
come from fermions with momentum perpendicular to 
q, and therefore might not be determined solely by the 
backscattering amplitude. 

The computation is in fact most straightforwardly car- 
ried out by rotating the s integral in Eq. 1(15(1 . leading in 
D = 2 to the non-analytic entropy per particle 



$Sna — 
3naCO = 



dT 



2vpkp 



dx xty(ix) 



(63) 



where x = v* F q/i\£ m \ and 
*(*) = 



dy $(y) 



(64) 



is the function whose discontinuity across the real axis is 
equal to <f> from Eq. 1 (15(1 . 

Nonanalytic contributions to Ea. H63|) can only arise 
from the dependence of the integral over x on its up- 
per limit, which is determined by the behavior of at 
large argument. Writing 



*(z — > 00) 



z 



yields 



Sna(T) 



* 5 



A(v* F k F y 



In 



>j f k f 



(65) 



(66) 



The frequency sum generates regular terms determined 
by the upper limit of the summation; these form an an- 
alytic expansion of S(T) in powers of T 2 . However, be- 
cause of the logarithm, the sum also produces a universal, 
i.e., independent of the upper limit, term of order T 3 



T y £ 2 ln ^kF 



-2C(3)T 3 



(67) 



where . . . stand for regular terms. This universal term 
can be obtained cither directly, by using Euler-Maclaurin 
summation formula, or by evaluating the frequency sum 
as a contour integral. Substituting (|67|l into (|66|l and 
differentiating over T, we obtain 



6Sna = *2 



T 



v* F k F 



+ 0{T 3 ). (68) 



Expansion in Particle-Hole Pairs 

First, we observe that because each factor of Pna 
brings an extra factor of z" 1 , the only contribution to 
^2 involves one real excited particle-hole pair and the 
leading correction to T, and is thus proportional to the 
square of the reducible interaction. Proceeding as in the 
previous section, we obtain the analogue of Eq. 120(1 

= - J d9 p V(ix; 6 p )T(ix; 9 P , 9 P ) (69) 

and the analogue of Eq. 1(21(1 

r = J d6 k T k {9 p -6 k )RcT(z,9 k )T k (e k ~6 p ). (70) 



Writing 



9 p -e k ) = J2 T LCos[L(6 p -0 k )], (71) 

L 
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using (cos# — ix)^ 1 = —i J dX exp[— iA(cos 9 — ix)] (for 
x > 0), expanding the exponential in terms of Bessel 
functions and using known integrals gives 



L 



r L r L , 



J , 2x 2 Vl + ar 



L+L' 



L-L 



(72) 
n 2 



Evidently, the large-argument limit of Eq. (|73|l is simply 



EC- 1 ) 



L+L' 



r L r L > 



L,L' 



(73) 

We see that the nonanalytic term in the entropy is de- 
termined by the backscattering amplitude. Restoring the 



spin and charge components of T, we then obtain for the 
entropy per particle 



SSna — 



3C(3) ( T ' (rt BS + 3it BS )+0(T 3 ). 

_ (74) 



2 V Vpkp 



This coincides with the result in Ref. 127 



Computation in real frequencies 

For the sake of completeness, we also demonstrate how 
to obtain Eq. I|68|) by performing a computation directly 
in real frequencies, without rotating the e integral in 
Ea. (|15|) . Using Eq. l|17fl an d the two-dimensional ana- 
logues of Eqs. 1 (20) 1 and (J2TJ, we obtain in 2D for the 
entropy per particle 



<5<Sna = 



T 



Vpkp 



dx- 



sinh^ 



x Jo 



Upfcjr / XT 



We know from Eq. (JOSJ) that SS NA oc T 2 . As T 2 is an 
overall factor in Eq. ((75)) , the rest of the integral can can 
be evaluated at T = upon which the upper limit of z 
integral is set to infinity. Integrations over x and over z 
then decouple, and evaluating the integral we obtain 



5S 



NA 



6C(3) 



T 



R, 



(76) 



where 



R = J dzz J dn k J dn p [r fe (n k -n p )] 2 

xlnxP (z,n k )Re-p (z,n p ) . (77) 



The subtlety, which we mentioned in the overview of this 
section, becomes clear if we assume momentarily that 
T k a constant. In this case, the angular integrations 
in Eq. l|77f) are performed independently, and R reduces 
to [r fe ] 2 / oo dzzImn(z)Ren(z). As Eas. (|61l62|l show, 
Imll (z) and Rell (z) are finite on the opposite sides of the 
boundary of the the particle- hole continuum (for \z\ > 1 
and \z\ < 1, respectively). One might then conclude that 
<5>Sna = 0. However, the confluent singularities in i?ell 
and Imll require a more careful analysis. Retaining the 



dzz J G?n k J drip Imll (z, n k ) Rell (z, n p ) [r fc (n k • n p )] ' 

r 



(75) 



infinitesimal iS to II at intermediate steps one finds 
/ dzzlmU(z) Ren(z) 



dzz 



Im- 



1 



z + l 
dxlm- 



i 



=Re- 



1+iS \/z-l+i5 



=Re- 



1 



\Jx + i yjx + i 
dy 7r 
_ 4' 



1 



(78) 



With this in mind, we proceed with the case of an 
arbitrary vertex T k (n k • n p ). As before, we expand r fc in 
harmonics f r,, see Eq. l(7T)l . Substituting this expansion 
into Ea. (|77fl . we obtain after simple algebra 



1? 



dz 



de v 



d9 u 



o Jo Jo 
xlmP (z,6» k ) Rc-P (z,8 p ) 

E 

L,L' 



X 

E 

L,L' 



T L T L , [cos (L + 9 k ) cos (L+9 P ) + (£ H 
TiXl'Il ,l'- 



(79) 
L-)}, 



where L± — L ± L'. As R turns out to depend only on 
whether L±, is odd or even we can just consider the L + 
term and double the result. 

Consider first contributions with odd L + . For these 
contributions, as we will see, one can neglect the i5 
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term in V, i.e., replace IniP (z, 6*k) by a 5— function 
Imp (z, 6*k) = — (l/2z)<5 — cos^k). Integrating then 
over z, we obtain 



h+L'=2P+l — — 
ir 



*' 2 cos[(2P+l)fl k ] 

«Pk 7 

o cos 6> k 



(80) 



The integral is convergent and equal to (— l) P 7r/2. Ac- 
cordingly, 



L+L'=2P+1 — 



(81) 



If L + L' is even, the integral over z vanishes unless one 
keeps 5 finite and takes into account that both Re P 
and Im P diverge as z approaches 1. This divergences 
come from the angular integral over a narrow range near 
9k = p = 0. Accordingly, we can safely set Ok = p 
outside the product Re V xlm V in Ea. (|79Jl . After simple 
algebra we then obtain 

7T 

Il+L'=2P = — T- ( 82 ) 



Combining Eqs.JSlJ and l|82|) . we obtain 



# = -j£(-i) i+L 'fLiv 

LX' 



TT 

4 



!? M)LfL ) 



(83) 

Substituting this result into Eq. (|76|l and introducing 
the charge and spin components of T, we reproduce Eq. 



Ring Diagram (RPA) Approximation in 2D 

In this subsection we present the two-dimensional ver- 
sion of the large-N calculation given above. Our goal here 
is to show that only Tg S determines the non-analyticity 
in S(T), whereas higher powers of Tl do not appear in 
the prefactor for the T 2 term in S(T). 

We again consider a simplified model with the effective 
interaction U = U +Ui cos(<^ p — fa), and spin degeneracy 
N = oo. The two dimensional version of Eq |44| is 



Sring — 



6C(3) 



T 



v* F k F 



J zfMmlnDet [1 - T(fa - fa)V{z, fa)} , (84) 



where now T = Ukp / (ttv f ) and 

1 cos fa 



V{z,fa) 



2tt cos ( 



(85) 



To evaluate the determinant we again must find the 
eigenvalues of the operator 1 — TV, i.e. solve 

AA(z, n p ) = A(z, tip)- J dn k T (n p ■ n k ) V(z, n fe )A(z, n k ) 

(86) 



We find the eigenfunctions of I— VP by making the ansatz 
A = A + Ai cos0 + A 2 sin0. We find that Det = £>i-Di 
with Di y2 given by Eas. l|46l47|l as before, but now instead 
of Eas. l|48l49|) we have 



Explicitly 



I'nm — 

hi = 

/oo = 1 
Ioi 

111 



(cost 



\n+m+l 



27T COS <j> — Z -1 

d<f) (sin0) 2 cos0 
2tt costb — z -1 



1 i 

z z\Jz 2 - 1 
1 1 i 



2 Z 2 Z 2^2— J 
- 1 1 Wz 2 - 1 

hl ~ 2'^ + ^ 2 ~- 
Evaluating and simplifying yields 

Di = (i + r)(i + ^) x 



• ffi i ri(i-r ; 

1 , r^i-rp) | *(r + -Ur-^ 



Vz 2 ^! 



with 



To 
fi 



i + r 



1 



Combining InDi and I11D2, we find 



(87a) 
(87b) 

(88a) 
(88b) 
(88c) 
(88d) 



(89a) 



D2 = I 1 + 4- I x I 1 - 4 + -^4 (89b) 



(90a) 
(90b) 



v' F k F /T rv* F k F /T 

zcMmlnDet[l - TV} = / zdzImlnQ(z) 



1 (fo + r, -?££!- H(i-r„) f^ + A 



VT- 



(91) 



The integral can be most straightforwardly evaluated by 
rotating the z— integral onto the imaginary axis z — » — iz. 
Then 



v* F k F /T riv F k F /T 

zdzlmlnQ(z) = — / zdzImlnQ(— iz) 

Jo 

(92) 
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As Q(—iz) is real and positive, the imaginary contri- 
bution from the integral over z can only come from 
the upper limit. One can verify that the dependence 
of the integral in i|92[) on the upper limit comes only 
from the 0(l/z) and 0(1/ z 2 ) terms in the expansion 
of ImhiQ(— iz) in powers of 1/z. Higher order terms 
in the expansion are ultraviolet convergent and can be 
safely evaluated by setting the upper limit to infinity. 
Expanding lnQ(— iz) to second order in 1/z, we obtain 

lnQ(- i z) = l-^^^^^ + 0(z^). (93) 

As in 3D, the leading term O (z^ 1 ) gives a renormal- 
ization of the effective mass, while the next term O (z~ 2 ) 
gives 




zdz lnQ(— iz) 



-(T - f 1 ) 2 Jmln[iv* F k F /T] = -(f - fx) 2 (94) 
Substituting Eq.J!2} into Eq.JS^J, we reproduce Eq. 



k 




k+q 



FIG. 4: A representative diagram for the thermodynamic po- 
tential, which contributes to the non-analyticity in the spin 
susceptibility to third order in the interaction. Vertical ar- 
rows denote spin projections for fermions in the corresponding 
bubbles. 



SPIN SUSCEPTIBILITY IN 3D 

In this section we discuss the relation between nonan- 
alyticities in the entropy and in the spin susceptibility, 
Xs- Until recently, the prevailing opinion had been that 
the nonanalytic T 3 In T dependence of the entropy is not 
paralleled by a nonanalyticity in Xs- Crucial evidence 
for this view was provided by the results of Carneiro 
and Pethick 0] and Beal-Monod et al. 0, who found 
that the leading term in the spin susceptibility scales 
as T 2 in 3D. However, in an important paper Belitz et 
al. |2lj demonstrated that the apparent analytic tem- 
perature dependence of Xs ma y be misleading. They 
performed a perturbative calculation of the momentum- 
dependent spin susceptibility Xs{Q,T = 0) at small Q 
and found a non-analytic Q 2 In Q behavior. Later, it was 
found [3(J that the magnetic-field dependence of a non- 
linear spin susceptibility parallels the Q-dependence, i.e., 
X s(Q = 0,T = 0,H)otH 2 \ii\H\. 

Nonanalyticity of the spin susceptibilit y w as also found 
for 2D systems by Millis and Chitov l24l and, later, 
Chubukov and Maslov ^2^, Galitski et al. |29j. and Be- 
touras et al These authors showed that Xs(T, Q, H) 
scales linearly with the largest out of the three param- 
eters (in proper units). Furthermore, Chubukov and 
Maslov [23 and Galitski et al. have shown that in 2D, 
the nonanalytic term in Xs may be expressed in terms of 
the spin component of the backscattering amplitude. 

These results call for a better understanding of the spin 
susceptibility of a 3D system. In the preceding sections, 
we found that the nonanalytic part of the specific heat 



consists of two contributions, the first one coming the ex- 
citation of a single particle-hole pair and the second one 
coming from the excitation of three particle-hole pairs, 
and with all coefficients determined by the harmonics of 
the Landau function. Here we ask if the non-analytic con- 
tributions to the spin susceptibility are given in the same 
way. It is clear from the previous second-order calcula- 
tions in 3D 

HJ0I22I 

that the single-pair contribution 
does contribute to the non-analyticity in Xs(Q, T, H), as 
is shown by, e.g., the second-order calculations. We show 
here that the three-pair terms, which were not considered 
in the previous literature, also contributes to the nonan- 
alyticity of Xs in 3D. For brevity, we refrain from a fully 
detailed analysis and simply indicate the origin of the 
effect. To this end, we consider the thermodynamic po- 
tential in a weak magnetic field perturbatively to third 
order in the interaction, which we take to be momentum- 
independent. We will be only interested in the spin effect 
of the field. A magnetic field (if) splits the Fermi surfaces 
for fermions with spins parallel and antiparallel to the 
field. This splitting does not affect the e/q nonanalytic- 
ity of the polarization bubble at small q, if a particle and 
a hole have the same spin (in this case, a magnetic field 
just shifts the chemical potential), but it has a nontrivial 
effect on a bubble composed of a particle and a hole of 
opposite spins, Hn(e,q;H). For a spin-invariant inter- 
action, such bubbles occur in diagram Fig. Q 3d. Label- 
ing the momenta as shown in the figure and integrating 
over fermionic momenta k, k', and k" and corresponding 
frequencies, we obtain the thermodynamic potential at 
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T = and finite H in the form 

S 3D {H) ~U 3 J de J d 3 qagne [ImU n (e, q; H)} 3 , (95) 
where 

Imn n {e, q;H) = - {-Ks/2q) 9 (v F q - \e-g L fi B H\) 

and and fis are the Lande- factor and the Bohr mag- 
neton, respectively. Note that Imll retains its H = 
form but the boundary of the continuum is shifted by 
the Zeeman energy. It is this shift which leads to a non- 
analyticity in 3(iJ). Integrating over q in Ea. (j95|l yields 

Z(H)~U S f de\s\ 3 \n\e-g L ti B H\. (96) 

Consequently, the nonanalytic part of S(-ff) behaves as 
H 4 In \H\ and thus Xs(H) cx H 2 In \H\. This contribution 
has the same functional form as the contribution from the 
spin component of the single-pair process. A similar con- 
sideration for a non-uniform field leads to a nonanalytic 
Q dependence of Xs(Q)- Xs{Q) cx Q 2 h\Q. 

Finally, at vanishing H and Q, both single-pair and 
three-pair processes contribute T 2 terms to the 3D spin 
susceptibility, which again is expressed only in terms of 
the Fermi-surface interaction and Fermi velocity. How- 
ever, it can hardly be identified in the experiment, as 
there are additional, purely analytic T 2 contributions to 
Xs{T) from all energies. 



is given by 

S = 4 T E / j^Z(urn,k)G (u; m ,k) (97) 

where Go(uj,k) — {iuj m — e^) -1 , and G _1 (ti;,fc) = 
G ! Q 1 (a;,fc) + E(w,fc). If the self-energy is independent 
of k (which is the case for, e.g., the electron-phonon 
interaction), the momentum integral involves only the 
Green's function, and J d D kGo(uj m ,k) ocsgnw. There- 
fore, 3 cx >0 E(aj m ). Evaluating the frequency 
sum, one obtains that the uj^ In \u> m \ behavior of the self- 
energy in 3D gives rise to a T 3 In T term in the entropy. 

In general, E(w m ,fc) depends on both uj m and k. We 
will see that in the low-energy limit, the self-energy in 3D 
is a sum of two terms, one depending on ilo — e&, another 
depending separately on ui m and 

E(w m , k) = (iw m - e fc ) + S (2) (u m ,e k ) . (98) 

On the Fermi surface (e* = 0), the two terms cannot be 
distinguished. However, away from the Fermi surface, 
their roles are different. If only the first term is present 
then the full Green function depends only on the variable 
uj — Ek so that the locations of the poles are not changed 
by the interaction and thermodynamics cannot change. 
Alternatively speaking, the first term in Ea. (|98|) affects 
only the quasiparticle renormalization factor Z, but as is 
well known, the Z factor does not affect the specific heat 
coefficient of a Fermi Liquid; hence the first term cannot 
contribute to the non-analyticity in S(T). 



FERMIONIC SELF-ENERGY 
Overview 

In this section we show how to obtain the non-analytic 
term in the entropy and the specific heat by evaluating 
the fermionic self-energy first, and then using the relation 
between the self-energy and thermodynamic potential. 
We present this calculation both because historically the 
nonanalytic terms in the specific heat were studied by 
first calculating the self energy and because, as we will 
show, the self energy exhibits several nonanalyticities, 
not all of which contribute to the specific heat. 

It has been known since the early days of the Fermi- 
liquid theory that at T — the subleading term in the 
3D fermionic self-energy on the Fermi surface E(w, fcj?) 
depends on the frequency in a nonanalytic way, as 
E(w, hp) oc u> 3 \nu>. Amit et al. 01 an d others argued 
that this oj 3 In u> dependence of E is the source for the 
T 3 lnT behavior of the entropy and the specific heat, as 
the entropy has the same scaling dimension as E and typ- 
ical frequencies lu are generally expected to be of order T. 
Indeed, to second order in the interaction U, the thermo- 
dynamic potential in terms of the Matsubara self-energy 



Explicit calculation of the fermionic self-energy in 

3D 



Now we proceed with an explicit evaluation of the non- 
analytic part of self-energy T,{u m ,k) in 3D. To simplify 
the presentation, we consider the self-energy only to sec- 
ond order in the interaction U(q), and only discuss small- 
angle scattering, i.e., assume U(Q) — U, L T (2fcj?) = [cf. 
Egs- flM - 

As we mentioned in Sec. II, the stationarity of S with 
respect to variations in G implies 0, ^| that one can 
neglect any explicit temperature dependence of E, i.e., 
it suffices to include only the T-dependence arising from 
the difference between T^2^ and J dui/2'K. 

The self-energy is given to leading order in U by 

k) = U 2 TY / J ^U(e m ,q)G ( 

(99) 

The computation is tedious but straightforward. The 
frequency-dependent part of the polarization operator for 
q <C kp is given by 



n(e m ,g) 



is,, 



2vpq 



■In 



It 



vpq 



IS,; 



VFq 



(100) 
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Expanding II (s m , q) in s m /vpq, we obtain 



of £( 2 )(w m . k) with the Green's function yields 



2v F q \v F qJ \ 



(101) 



Substituting Eq. i|l(Jl|l into Ea. 1991) . integrating over q, 
and summing over e m (which, to logarithmic accuracy, 
is equivalent to just integrating over s m ), we find that 
E(w m , k) is given by Eq. [HB] with 



£ (1) (^ m - efc) = 



mkF u) ilUJ ~j 3 In- 
384£ 2 



(102) 

This contribution originates from the first, leading term 
in the expansion of II (e m , q) in e m in (|101J) . and it comes 
from internal bosonic frequencies that exceed the exter- 
nal one, i.e., from |e m | > |w m |. The second term depends 
on uj m and ek separately and has contributions from both 
the first and second terms in Eq. (lOlj . To logarithmic 
accuracy, 



v(2)/ us ■ ( mUk F \ 

S v > {u m , k) = -l I — I 



2 ,,3 



16£| 



In- 



E F 



-O 



{iuj m - e fc )ln 



lW m + efc 



iw + e fc 

(103) 



(see comment .36]). Note that T,^(uj m ,k) comes from 
the range of bosonic frequencies between — |w m | and \ui m \, 
i.e., relevant internal frequencies which are smaller than 
the external one. 

Comparing Eas. ((102|l and 1)103(1 . we see that for a 
generic ratio of efc/w m and, in particular, for a parti- 
cle on the Fermi surface, where efc = 0, both contribu- 
tions behave as Sfc^m, kp) oc wj^ln \ui m \; this result is in 
agreement with Ref.|l2j. However, only £( 2 )(o; m ,fc) ac- 
tually gives rise to the T 4 lnT term in 2, and hence to 
the T 3 lnT terms in the entropy and the specific heat. 
Indeed, the contribution of EP'(w m , k) to 5 is 



T£ / de k E« (u! m ,k)G (u! m ,k) 
de k {iuj m - e fe ) 2 ln 



(iu m - e fc ) 2 



= 2TVRe / dCC 2 ln-f, C = e*-*Wm- 

(104) 

The resulting integral is determined by the upper cut- 
off of the integration, and is an analytic function of w m , 
hence gives only an analytic contribution to the thermo- 
dynamic potential. On the other hand, the convolution 



J de k JlW{w m ,k)G {w,k) 



oc 



de k 



hi- 



6k + 



(105) 



To logarithmic accuracy, the integral over e k gives 



de k ^ E F 



de k . 

In ■ 



= «7TSgnCJ„ 



Ef de k E F 
1 1 

\u m \ e k |W m | 



= 2z7rsgno; m In 



Ep 
\w m \ 
(106) 



Substituting this into Eq. I|97|) . evaluating the frequency 
sum using 



In 



E F 



= 2^ T 4 ln ^ (107) 
15 T V ; 



and differentiating over T, we obtain for the entropy per 
particle 



S(T) = -— 



lUkr 



T 



v* F kp 



In 



'>pK F 

T 



(108) 



This coincides with Eq. (|39|l expanded to second order, 
and evaluated at U(0) — U, U(2kp) = [in this ap- 
proximation, K c = mkpU/ir 2 , and K s = 0]. We also 
verified that the contributions to £( 2 ) from e m /q and 
e^n/q 2 terms in IT yield equal contributions to 5. This is 
consistent with our earlier observation that non-analytic 
term in 3 comes from the cross-product of the 0{e m ) and 
O(e^) terms in IT. 

We see therefore that non-analytic term in S(T) can 
indeed be obtained by evaluating the self-energy first 
and then convoluting it with the fermionic Green's func- 
tion. However, only a part of the self-energy actually 
contributes to the nonanalyticity in S(T). The term in 
the self-energy that depends on iui m — e k gives a renor- 
malization of the quasiparticle weight and does not lead 
to a nonanalyticity in S (T) . 

It is instructive to comment here on the role of for- 
ward scattering. In 2D, forward scattering is special in 
the sense that its contribution to the self-energy is com- 
parable to that from backscattering [53; both are of 
order e 2 ^ In |e m |. In fact, the 2D analogue of £W(e m , k) 
comes entirely from forward scattering. In 3D, forward 
scattering is much less effective because of phase space 
restrictions. One can show that in this case T,^(e m , k), 
evaluated to all orders in the perturbation theory in the 
dimensionless coupling u, is not restricted to forward 
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scattering but rather comes from a wide range of scatter- 
ing angles. As we will discuss in the next Section, there 
still exists w 3 logw contribution to H(uj 7 k F ) specific to 
forward scattering in 3D, but this contribution is expo- 
nentially small for a weak interaction. 



COMPARISON TO EARLIER STUDIES 

As we said in the Introduction, the analysis of the 
T 3 In T term in the specific heat of a 3D Fermi Liquid 
has a long history. In this section we compare our find- 
ings with the existing literature on the subject. We be- 
gin with reviewing the T 3 In T nonanalyticities for the 
cases of electron-phonon and fermion-paramagnon inter- 
actions, and then proceed with the generic Fermi-liquid 
case. 



Electron-phonon interaction 

The electron-phonon coupling leads to nonanalytici- 
ties in the entropy; in fact, the nonanalytic correction 
to the Fermi-liquid fixed point was discovered first by 
Eliashberg || in a study of acoustic phonons coupled to 
fcrmions via a deformation potential. Integrating out the 
phonons, one obtains a new contribution to the electron- 
electron interaction. For optical phonons, or for acoustic 
phonons with a piezoelectric coupling, this new interac- 
tion is frequency dependent on the scale of the Debye 
frequency, but at low frequencies it just generates an ad- 
ditional contribution to the electron vertex r fc , so the 
results we have presented in previous sections carry over 
directly. 

The case of an acoustic phonon-deformation potential, 
studied by Eliashberg, presents a new feature: the ef- 
fective interaction itself contains a nontrivial e 2 /q 2 term 
which leads to a new non-analytic contribution to entropy 
in D = 3 (but not in D = 2). We discuss this case in 
more detail. Perturbation theory in the electron-phonon 
interaction is controlled by a small Migdal parameter - 
the ratio of the sound velocity c to the Fermi velocity vp. 
To leading order in c 2 /y 2 ?. the thermodynamic potential 
can be expressed as 0, l37| 



^e—ph 



\ J ^n B (e/T) J Aim [mi?- 1 (e, q)] 



(109) 

where D 1 is the inverse phonon propagator given by 

D- 1 (e, q) = Dq 1 (e, q) + 2g%c? ] q 2 U(e , q). (110) 

with bare propagator 

D 1 (E,q) = -E 2 + c 2 q 2 . (Ill) 

Here <7 p h is the dimensionless effective electron-phonon 
coupling. In a typical metal, g p h ~ 1 0. The expres- 
sion for the thermodynamic potential 3 e _ p /j is similar 



to the RPA expression for electron-electron interaction, 
c.f. Eq. (|42|l . but there are two importand differences. 
Firs, the bare propagator D describes gapless excita- 
tions. Second, the second term of Ea. (|110l) contains an 
extra factor q 2 . This factor guarantees that the q = 
phonon mode does not affect the electrons ( "Adler prin- 
ciple"). These two differences combine to preserve the 
T 2 logT form of the nonanalytic correction to 7 in 3D. 
One can rearrange the Im In . . . in Eq. Ijl09|l by factoring 
out purely real quantities to get 

=.2 



Im In D 1 = Im In 



1 - 



)V 



S*n(e,g) 



(112) 



where now (c*) 2 = c 2 + 2g 2 h ReIl(q, 0) is the sound veloc- 
ity renormalized by the analytic terms in II, and g* sim- 
ilarly is the renormalized coupling, g 2 = 2g 2 h (c/c*) 2 . In 
analogy to the previous section, we now expand in powers 
of Im!!. As before, the term with one real particle-hole 
pair (i.e., of first order in Imll) combines with the lead- 
ing e 2 /g 2 term in the interaction vertex to give a T 3 InT 
nonanalyticity. In addition, however, another T 3 logT 
emerges from the combination of Imll and ~ e 2 /c* 2 q 2 
term from the bare phonon propagator. Such term does 
not exist for electron-electron interaction. A simple anal- 
ysis shows that this new contribution is larger by the 
inverse Migdal parameter (vf/c) 2 than the contribution 
from the interaction vertex and hence is the dominant 
phonon T 3 log T contribution to the specific heat 

Note that the sign of the phonon T 3 In T term is op- 
posite to that of the electron one, i.e., d~f(T)/dT is pos- 
itive for electron-phonon contribution and negative for 
the electron-electron contribution 

In D = 2 (when both electrons and phonons are two- 
dimensional) the q dependence of the phonon propagator 
remains the same, but the integration measure changes, 
and g 2 h term does not give rise to a nonanalyticity in 
S. Only the "conventional" term of order g^ h gives rise 
to non-analyticity. From this perspective, the electron- 
phonon problem in 2D does not differ qualitatively from 
the electron-electron interaction. In the case of mixed 
dimensionality, e.g., planes of electrons embedded into a 
3D elastic continuum, more complex situations are pos- 
sible. 

One may also extract the "phonon" T 3 In T term in the 
entropy from the nonanalytic form of the electron self- 
energy for the electron-phonon interaction Eph(w m ) ~ 
o^ln |u; m | 0|, in analogy with Ea. (|9"?|) . We re-iterate 
that a direct relation between the nonanalyticities in the 
self-energy and entropy exists only for the case when the 
self-energy depends only on frequency. This point will 
be important for the discussion of the zero-sound mode 
contribution to the entropy in Sec B 

Finally, we note that the large phonon-induced non- 
analyticity in the entropy is not paralleled by a similar 
term in the spin susceptibility. This comes about be- 
cause phonons contribute directly to the charge vertex, 
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and only inderectly to the spin vertex, at the subleading 
order in the Migdal parameter. The basic reasonin g fo r 
this was given long time ago by Prange and Kadanoff [3!j , 
who showed that in the Migdal approximation, phonon 
contribution to the thermodynamic potential is confined 
to the Fermi surface, and does not change appreciably 
when the Fermi surface position changes. As a result, 
an applied magnetic field, which shifts spin up and spin 
down Fermi surfaces, does not change the phonon contri- 
bution to the thermodynamic potential, to leading order 
in the Migdal approximation. Mathematically, this can 
be seen as follows: (i) The static spin suscepbility Xsil) is 
given by the particle-hole bubble, (ii) For free fermions, 
Xs(q) comes from the states away from Fermi surface, 
if one integrates first over fermionic dispersion e^, while 
the contribution from the states near the Fermi surface 
vanishes because the poles of the integrand as a function 
of dispersion ek reside in the same half-plane, (iii) The 
renormalization of the spin susceptibility due to electron- 
phonon interaction occurs via self-energy insertions into 
the spin-polarization bubble. As the self-energy does not 
depend on the electron momentum, these insertions still 
leave the two poles of the integral over in the same 
half-plane, hence the integral over low-energies still van- 
ishes. 

The analysis beyond leading order in the Migdal 
parameter requires further analysis, not given here. 
However, for comparable electron-phonon and electron- 
electron couplings, the phonon renormalization of Xs is 
of the same order as the purely electron one. 

Paramagnon model 

In the paramagnon model 10], the fermion-fermion 
interaction is approximated by an interaction between 
fermions and overdamped long-wavelength spin fluctu- 
ations (paramagnons) The model has been extensively 
used to describe itinerant electrons near a ferromagnetic 
instability. While the analysis is almost identical to the 
analyses we have already presented, the possibility tun- 
ing the model through a ferromagnetic quantum phase 
transition raises a new issue which requires discussion. 

The RPA thermodynamic potential per unit volume in 
the paramagnon theory has the same form as Eq. (|42(l but 
with slightly different normalization factors; we write it 
explicitly here for convenience: 

E = Ea + l J V ns(£/T) / (£p (113) 

Here So is the thermodynamic potential of the nonin- 
teracting fermions and Dg 1 (e,q) = 1 + 2gV p h(£,q) is 
the inverse propagator of spin fluctuations, g is the spin- 
fermion coupling, and V p h is the full polarization bub- 
ble for noninteracting fermions. This model is therefore 



equivalent to the previously studied RPA case, but with 
only the L — spin channel interaction retained, so that 
the determinant [cf. Eq. 144fl ] may be written down im- 
mediately. However, one extra piece of physics becomes 
important. This theory can describe a ferromagnetic in- 
stability, and indeed an original motivation for the model 
was the study of effects of long wavelength spin fluctua- 
tions, which may be expected to be important in nearly 
ferromagnetic materials. Near a ferromagnetic instabil- 
ity, the dependence of the analytic part of n on q becomes 
important. Decomposing n into the analytic (const + q 2 ) 
and nonanalytic (Hna) parts, we have 

D^faq) = (l + B)+e o q 2 + 2gIi NA {e,q) (114) 

with B = 2gll(0, 0) — > —1 as the ferromagnetic transi- 
tion is approached. We now proceed as before expand- 
ing in lmLljvA (representing real particle- hole pairs) and 
RenNA (representing the nonanalytic part of the dynam- 
ical interaction. The expansion contains inverse powers 
of 1 + B + £o<Z 2 : so that the power counting of the mo- 
mentum integrals changes for q > £ = £o/\A + B |4C| - 
The nonanalyticities we have discussed in this paper arise 
only from the Fermi- liquid regime T/vp < q < £ . As 
criticality is approached, £ diverges and the temperature 
window in which the T 3 lnT term may be observed be- 
comes vanishingly small. It is thus somewhat misleading 
to state that the T 3 In T nonanalyticities become the crit- 
ical nonanalyticities; rather, they should be regarded as 
a property of the Fermi-liquid regime only, which is the 
"quantum disordered regime" of the ferromagnetic quan- 
tum phase transition. 

With this proviso, the previously given analysis can 
be applied directly. We note here that the original work 
0,[ljj focussed on the three-pair contribution, obtained 
by treating ReHjVyi as a constant. In 3D one finds 

(!)' G£)'-^ (-» 

with B = -235(0,0) and A = (£0/£) 2 = 1 + B. This 
expression coincides with Eq. I|39(l , if one considers only 
spin contribution, neglects the backscattering term in K s 
and all F except for r Sj o, and identify B with the Landau 
parameter T Sy o, such that F Sj o = B/(l + B). Note that 
the narrowness of the "Fermi-liquid" momentum regime 
appears here as a dependence of the upper cutoff of the 
logarithm. Incorporating the cutoff into the analytical 
T 3 terms, while mathematically justified, obscures this 
physics. 

The single-pair contribution, overlooked in the param- 
agnon literature |lfl Il7| , arises (as previously discussed) 
by including the nonanalytic momentum dependence of 
Reli and is 

(!)'G£) ,ta ^- <»«> 
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This expression also coincides with Eq. (|39[1 , if one again 
considers only the spin contribution, but now takes into 
account only only backscattering term in K s , neglect all 
r except for T s q 1 and again identify B with the Landau 
parameter r Sj o- Comparing the two non-analytic con- 
tributions to 6S(T), we find that for negative B, i.e., 
positive g, they have the same sign and become equal 
at B rs —0.55. Further away from the ferromagnetic in- 
stability (at smaller 1 + B), the single-pair contribution 
dominates, while closer to the instability, the three-pair 
contribution is larger. 

The same situation holds in 2D. Away from critical- 
ity, there is only one nonanalytic contribution to S(T), 
from H 2 (e m ,q), i.e, from a single-pair excitation with a 
frequency-dependent interaction. It yields 



T 



(117) 



Eq |117| only applies for T < Vp/£; again as criticality is 
approached the window of validity closes. 



Fermi Liquid 

The nonanalytic term in the entropy of a 3D Fermi 
Liquid arising from a generic fermion-fermion interaction 
was considered in a number of publications in 60s and 
70s 0, Q, S [H Q El E3 • Recently, the T 3 In T term 
in S (T) was reproduced via multidimensional bosoniza- 
tion ME!. 



Zero-sound mode 

The very early studies d, 0, Q focused on the non- 
analytic contribution due to the interaction with zero 
sound. As this theme returns in later work (see, e.g., 
Ref . [l^| ) . it is worthwhile to summarize the status of the 
problem here. The idea that the interaction with zero 
sound may lead to a nonanalyticity in the specific heat 
was put forward by Anderson || . Following this sugges- 
tion, Balian and Fredkin Q considered a phenomenolog- 
ical model which treated the interaction with zero sound 



J 



in analogy with the electron-phonon interaction. How- 
ever, an unphysical choice of the interaction vertex (it 
remained nonzero at q — > 0, and thus did not satisfy the 
Adler principle) led Balian and Fredkin to the conclu- 
sion that the nonanalyticity in the self-energy was ver 
strong (w m In |w m |). Later, Engelsberg and Platzman 
pointed out that the zero-sound vertex vanishes as q for 
q — ► 0, as it does for the deformation-potential- type inter- 
action with phonons [cf. Eas. (|109lll0(l ]. They obtained 
an u;^jln|a; m | non-analyticity in the self-energy at the 
Fermi surface, with prefactor proportional to the differ- 
ence of the zero-sound and Fermi velocities c zs — vp > 0. 
In 3D, this difference is exponentially small for weak in 
teractions u: c zs — vp oc e 
analysis yields an ui^ In \u m \ 
a prefactor of order it 2 . 



l / u . In 2D, the analogous 
term in the self-energy, with 



In this and all subsequent publications, which em- 
ployed the zero-sound nonanalyticity in the self-energy, 
it was assumed that this automatically implies a nonan- 
alytic entropy. However, the two of us and collabora- 
tors have shown recently [2?l Eij ] that this assumption 
is, generally speaking, incorrect, and that the interaction 
with zero sound in 2D does not lead to a nonanalytic- 
ity in the entropy to any order in the interaction. The 
most general argument is that the thermodynamic poten- 
tial can be obtained directly in Matsubara frequencies, 
by expanding the Luttinger-Ward functional in the bare 
(fermion-fermion) interaction. The zero-sound mode sim- 
ply does not arise in this formalism. Recently, Catelani 
and Aleiner [42j argued away the forward-scattering (and 
thus the zero-sound) contribution to thermodynamics on 
the basis of gauge invariance of the Fermi-liquid kinetic 
equation. 

On a more technical level, we come back to the argu- 
ment that a direct relation between the nonanalyticities 
in the self-energy and entropy exists only if fermions in- 
teract with slow boson modes, so that the self energy 
is approximately fc-independent ("local"). As the zero- 
sound velocity is necessarily larger than the Fermi veloc- 
ity, this relation does not have to hold for this case. The 
breakdown of this relation becomes evident on examina- 
tion of the real- frequency analogue of Eq. (|§7j) 



SS = 



2 d 




r d D k 


poo 


TT dT 


T j 


1 (2tt)D 





duiLo^^ {Re£(w, k)TmG (u, k) + Im£(w, fc)ReG (w, k)} 



(118) 



r 



where Go(u>, k) — (uj — ek + i0 + ) and np is the Fermi func- 
tion. In a local theory, the k integral of the ReSImGo 
term in Ea. (|118|l vanishes, and only ReS contributes to 
SS. In a non-local theory, both the ReS and the ImS 



terms contribute. An explicit computation shows (27l l28| 
that the zero-sound contributions from the two terms in 
Ea. pi8[l cancel each other in 2D. Although we have not 
performed such a calculation in 3D, we believe that the 
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result will be the same as in 2D, i.e., the zero-sound con- 
tribution vanishes. This also follows from the simple ar- 
gument that in the RPA (ring) approximation in 3D, the 
entropy for the interaction with zero-sound is given by 
Ea. (|42[l . in which the argument of the logarithm is re- 
placed by the inverse effective interaction U e g(e, q). Near 
the zero-sound pole, U e g(e, q) oc e~ x l u q 2 j [e 2 — c 2 s q 2 ) 
and the zero-sound contribution to the entropy coincides 
with that of a free bosonic mode. Accordingly, the en- 
tropy is analytic and scales as T 3 0, ^| . 

Fermion-Fermion interaction 

Detailed calculations of the self-energy and S(T) in 
generic models of interacting fermions have been per- 
formed by Amit, Kane and Wagner an d Pethick and 
Carneiro ^J, and we compare our results with theirs. 

Amit and co-workers computed the self-energy near 
the mass shell, expressed the result in terms of Lan- 
dau parameters keeping only the partial amplitudes with 
L = 0, 1, and then used the resulting expression for the 
self-energy to compute the T 3 In T correction to the en- 
tropy obtaining a result somewhat similar to the one pre- 
sented here. There are however three qualitative distinc- 
tions between their results and ours. First, Amit et al. 
argued that the non-analytic In \ujm\ term in the self- 
energy comes only from the Landau damping term (e/q) 
in the polarization bubble. As we have shown, there are 
actually two contributions, one from the Landau damp- 
ing (e m /q) term in II( m , q) (considered by Amit and co- 
workers) and another from E 2 m jq 2 (overlooked by Amit). 
Second, Amit et al. argued that the non-analytic self- 
energy at the mass shell partly comes from the interac- 
tion with zero sound. We showed in the previous subsec- 
tion that this interaction does not contribute to T 3 In T 
term in S(T). Third, Amit et al. argued that 2kF scat- 
tering does not contribute to thermodynamic potential 
to second order in the interaction. We found that the 
contributions from small momentum scattering and from 
2kp scattering are identical at this order (up to prefac- 
tors U (0) and C/(2fc,p)). We also note that the final result 
of Amit et. al. for 5S(T) [Eqs. (IV.17) from 0b and 
(VI. 8) and (VII. 18) from |l2ja is not identical to the re- 
sult we presented. Their combination of Aq and A\ does 
not reduce to the square of the backscattering amplitude, 
and their overall factor in SS(T) and the relative prefac- 
tor for the terms with 7r 2 /12 in K s in (|39[) are different 
from ours. 

Our result for the entropy fully agrees with that of 
Carneiro and Pethick if we use the same ap- 

proximation as they did, namely, neglecting all Lan- 
dau coefficients with L > 2. Within this approxi- 



mation, our Eq. l|39|l coincides with Eqs. (22) and 
(A19) in Ref.[l5(. In particular, to second order in the 
scattering amplitude, their prefactor for the entropy is 
Al + A\ + A\ + 2A A 2 - 2A 2 A X - 2A Ax (the last term 
is missing in (A19) of Ref.|15|. but this is obviously a 
misprint). This combination is nothing but the square 
of the backscattering amplitude ^4| s = (A — A\ + A 2 ) 2 , 
which is the first term in our prefactor K in l)39|) . From 
this perspective, our new result is the observation that 
sum of all bilinear products of partial components of A 
reduces to the square of the backscattering amplitude. 

It is frequently stated in the literature that Pethick and 
Carneiro considered forward scattering between Landau 
quasiparticles, apparently in disagreement with our re- 
sult that forward scattering does not give rise to a non- 
analyticity in 5. There is in fact no contradiction be- 
cause the forward scattering considered by Pethick and 
Carneiro involves small angle scattering between parti- 
cles slightly away from the Fermi surface, whereas we 
consider scattering of particles at the fermi surface. 

To see this explicitly we note that Pethick and Carneiro 
used the Fermi liquid relation between the entropy and 
the thermal correction to the quasiparticle energy 

5S(T) = "£Ae p (T)-^, (119) 
p 

where n p is Fermi function (a sum over spins is implicit), 
expressed Ae p (T) in terms of the Landau function 

A6 p (T) = ]T/ p , p+q n p+q , (120) 

and argued that when / is expanded in powers of p-q, 
the expansion contains a quadratic term (p-q) 2 . This 
quadratic term gives rise to a T 3 In T contribution to 
AS(T). 

It is important to realize that the Landau function 
/p,p+q describes the interaction between physical quasi- 
particles, for which a; = v* f (\p\—pf)- Accordingly, / p , p + q 
coincides, up to an overall factor, with the fully renormal- 
ized vertex r[(p,w), (p + q, w + e); (p'w'), (p' + q, u' + e)] 
(where the first two pairs of arguments corresponds to the 
initial states and the last two to the final ones), taken at 
p = p',u' = uj and e = e p + q — e p ~ vftci p ■ q, where 
n p is a unit vector along p. In fact, the (p-q) 2 term 
in /p,p+q describes how the interaction between quasi- 
particles evolves when they move away from the Fermi 
surface. A similar consideration has recently been ap- 
plied to the analysis of nonanalytic terms in the spin 
susceptibility pjf. 

The full irreducible vertex 
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r[(p,w), (p + q,uj + v F n p • q); (p'w'), (p' + q, J + v F n p ■ q)]| e=UF „ p . q = T Prpl (9 ), (121) 



where 9o is the angle between p and q, satisfies an inte- 
gral equation, graphically shown in Fig. 1 of [l5|. The 
quadratic dependence on cos#o, which eventually leads to 
T 3 lnT term in S(T), comes from a virtual particle-hole 
bubble, composed from quasiparticles at finite energy, 
i.e., away from the Fermi surface. All analytic correc- 
tions are absorbed into the bare vertex r pp , (cos^o); m 
which one can then safely set cos#o = 0. Explicitly, 

r p , P '(0o) =r° )J) / +A J d 3 pir p Pi G Pl Gp 1+q r Pli p/(<7, 6» ), 

(122) 

where A is a constant prefactor. The bare vertex de- 
scribes the interaction between physical fermions at the 
Fermi surface and coincides, up to an overall factor, with 
r fe (n p • n p /) from Eq. 0- However, r PiP (6') cx / p , p + q 
is only obtained by solving the integral equation (|122Jl . 

Pethick and Carneiro solved Eq. (|122Jl keeping only 
L = 0..2 partial components of T k and found contribu- 
tions to cos 2 9o term in r p p (6'o) of the second and third 
orders in T p p , . The second-order contribution and a part 
of the third-order one correspond physically to processes 
with one excited particle-hole pair, while the rest of the 
third-order contribution describes the process with three 
excited particle-hole pairs. 

One can easily demonstrate, starting from Eq. i|122fl . 
that the second-order contribution is indeed expressed 
in terms of the backscattering amplitude, as we found 
in Eq. (|39|1 . Applying one iteration to Eq. Ijl22|l ancl 
integrating over e Pl and over frequency, we obtain 



F 



p,p(V ) cx cos6> 



LIZ LQ P,Pl Pl,P 



(123) 



47r cos 8q — cos 9 

where fig is a solid angle and 9 is the angle between p! 
and q. Substituting an expansion of T p =T°(9 — 6>o) 

(124) 



cos 



-o ), 



COMPARISON TO EXPERIMENT 

In this section we present a brief comparison of the re- 
sults to measurements of the T 2 In T term in the specific 
heat coefficient 7 (T) = C(T)IT in 3 He Q and in sev- 
eral heavy-fermion materials 0,^3- The 3 He data were 
analyzed using a formula equivalent to our Eq. I39l but re- 
taining only L = 0, 1 harmonics in charge and spin chan- 
nels. The values of r c .o, L Cj i, and T s ^ were taken from 
independent measurements of the compressibility, effec- 
tive mass, and spin susceptibility, respectively, and r Si i 
was chosen as a free parameter to fit the observed T 2 In T 
term in7(T). At zero pressure, r Cj0 = 9.15,r Cj i = 5.27, 
r s , = -0.700, and T sA = -0.55 @|. Use of Eq.O 
then gives Tbs.c = —1-01 and Tbs,s = —1.66, which 
yields for the combination of backscattering amplitudes 
entering Eq.^ L2 BS ,c + 3L2 BS:S = 9.28. Estimating 
the rest of the terms in K c + 3K S (which are cubic in 
the amplitudes) in the same approximation, we obtain 
(K c + 3if s ) cubic = 20.05. Thus in 3 He at ambient pres- 
sure it appears that the backscattering contribution is 
about a third of the total result. At higher pressures, 
the relative importance of the backscattering contribu- 
tion diminishes to about a quarter of the total result. 
This suggests that in 3 He the "paramagnon" (r 3 ) con- 
tribution is quite important, although the lack of infor- 
mation about higher Landau parameters means that this 
conclusion must be regarded as provisional. 

Similar analysis of the data for UPt3 was performed in 
Ref. ^ij- In this case, however, only the L = harmon- 
ics in Ea. l|39|l were retained. The value of r Sj o extracted 
from the magnitude of the T 2 InT term in 7(T) does not 
agree with the value obtained from the enhancement of 
the spin susceptibility. Given this controversy, we refrain 
from estimating the relative importance of the backscat- 
tering contribution to j(T) for this system. 



into Eq. (|123(l , expanding the integrand to order cos 2 9o 
and evaluating the prefactors in the same way as in Sec. 
IV, we obtain after rather straightforward algebra that 



r PiP (0 o ) oc cos 2 9 Q f „f m (-l)"+ m = cos 2 8o T 



2 

BS 



(125) 

This fully agrees with Eq. (|3*9"j> . 

The same analysis can also be applied to the 2D case. 
The only difference is that in 2D the T 2 term in S(T) 
comes from all angles 9$, rather than from a specific range 
of small 9q. Still, one can obtain an integral equation for 
r PiP /(6'o), similar to Ea. i|122|l and solve it by the same 
method as in Sec. V. 



SUMMARY 

We have revisited the old subject of nonanalytic con- 
tributions to the entropy (or specific heat coefficient) in 
light of recent progress in our understanding the struc- 
ture of the leading corrections to Fermi- liquid theory. We 
have clarified the connection to the Fermi Liquid fixed 
point by showing that the coefficients of the nonana- 
lytic terms may be expressed in terms of partial har- 
monics of the fully reducible static scattering amplitude 
r fc (n p ,n p /) with both momenta taken on the Fermi sur- 
face. We have presented a formalism which allows us to 
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distinguish between excited particle-hole pairs and their 
(possibly nonanalytic) interaction. 

We have found two classes of contributions to the spe- 
cific heat nonanalyticity. The first one, pertinent to both 
two- and three-dimensional Fermi liquids, arises from an 
excitation of a single-particle pair above the ground state 
combined with a nonanalyticity in the quasiparticlc inter- 
action vertex. The second one, pertinent only to Fermi 
liquids in dimension D > 3, arises from the excitation 
of three particle-hole pairs, which interact via the ana- 
lytic fixed-point interaction. In D < 3 the first sort of 
nonanalyticity is determined solely by the backscatter- 
ing amplitude, but in D > 3 additional contributions 
may arise. 

The fact that the non-analyticity in the entropy in- 
volves the spin component of the scattering amplitude 
suggests that the same physics is responsible for the non- 
analyticities in the susceptibility which have been dis- 
cussed extensively elsewhere 0, HjJ H3, H3| • We have 
presented simple calculations which confirm this sugges- 
tion. We thus conclude that in D = 3 the nonanalytic 
momentum and magnetic field dependences of Xs, re- 
ported earlier, involve terms up to third order in the ex- 
act quasiparticle scattering amplitude, and the quadratic 
term is the square of the spin component of the backscat- 
tering amplitude. From this perspective, the nonanalyt- 
icities in the spin susceptibility are the same nonanalyt- 
icities which have been studied in the context of specific 



heat for many years. 
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APPENDIX A: 2k F CONTRIBUTION TO 
NONANALYTICITY IN THE SPECIFIC HEAT. 

We argued in Sec. II that one can can re-express the 
2kp contribution to the thermodynamic potential 5 as 
an effectively small q contribution. For completeness, 
however, we show in this Appendix how to evaluate the 
2kp contribution to 3 directly. The technical difficulty 
of this approach is that one needs to know the polar- 
ization bubble near 2kp at finite frequency and at finite 
temperature. For simplicity, we limit our analysis to the 
second-order perturbation theory. 

We first obtain the asymptotic form of II (e m) q, T) and 
then use it in the calculation of the thermodynamic po- 
tential. It is convenient to separate H(e m ,q,T) into a 
T = part and the T— dependent part as 

n{e m ,q,T) = n + n T , 

where 



/d 3 k 
(2tt ) 



/d k 
—3 [n F (e k ) - 6 (-e k )] 
(2tt) 



d 3 k 
(2nf 
d 3 k 
(2ir) 3 



1 



1 



lS m — £k+q + £k l£m — e k + £k-q. 
1 



lE m — £k+q + £k Km — £k + e k-q. 



and hf (ek) = 1/ (exp (e k /T) + 1) is the Fermi function. 
Consider the T = part first. Expanding fermionic dis- 
persions near q « 2k p yields 



£k±c 



-e k + v F q + 2E F a 2 



where q — q — 2kp <C 2kp and a <C 1 is the an- 
gle between — q and k (for ek+ q ) and between q and 
k (for 6k- q ). We assume that the k integral is con- 
fined to the Fermi surface, and replace J d 3 k/ (2ir) 3 • ■ • — > 
(vp 12) J dek J daa . . . , where vp = mkp /2tt 2 is the den- 
sity of states at the Fermi level per one spin orientation. 
Integrating over e&, we arrive at 



Vp 

n = — / ax In 



(0 



n 2 



where 

f2 m = e m /Ep,Q = vpq/Ep, and x = 2a 2 . 
Performing an elementary integration, we obtain 

Up 



n = — 

16 



7r|fL 



(126) 



2|f2 m | arctan 



Q_ 



QHn 2 m + Q 2 ) 



Notice that since we restricted the momentum integra- 
tion to the vicinity of the Fermi surface, our expression 
does not contain a regular part of II, which comes from 
all energies. However, in what follows we will need only 
a nonanalytic part of II. 

At e m — 0, Ea. H126fl reproduces the static Kohn 
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anomaly 



n 



0|e m=0 



VF 

' 8 



Qln|Q|. 



(127) 



(128) 



At finite e m and q = 2Uf, i.e., Q = 0, 

vp - 

no| g =2fc F = — -j^rl^ml- 

This expression implies that the spectral weight of 
particle-hole pairs with the center-of-mass momentum 
near 2k f is proportional to the excitation frequency. 

For the T— dependent part, we find after integrating 
over a, 

n T = -f r^n F (^ F )ln ^ + ^ + % , (129) 

where £ = e^/Ep- 

The 2k F contribution to the thermodynamic potential, 
5 2 fc F , is given by 



E 2kF = -ZNU'T 



(130) 



^JdQ [n 2 + 2n n T + n|] 



where N is the number density. The choice of the cutoff 
in the Q integral does not affect the result to logarithmic 
accuracy. The rest of the calculation is simplified by 
making two observations. First, the first term (|f2 m |) 
in Eq. I|126|> is independent of Q, while the rest of the 
terms in this equation as well as II t are odd in Q. The 
square of the first term in Eq. (|126[) does not lead to a 
nonanalyticity in S, whereas the cross-products of this 
term with the rest of II is odd in Q and thus vanish upon 
integration. Therefore, the first term can then be safely 
eliminated from Eq. (|126f) which can then be written as 



16 

2|O m I arctan 



Q_ 



QHn 2 rn + Q 2 ) 



(131) 



Second, typical £ in Eg . ((H^ are of order T = T/E F , 
while the In T dependence of the thermodynamic poten- 
tial comes from the region Q 2 ,f2^ ^> T 2 . The integral 
in Ea. p29|l can then be expanded in powers of £, and a 
simple analysis shows that one needs to keep only terms 
of order £ and of order £ 3 . With this simplification, 11^ 
becomes 



tt 2 vf 
12 



T' 



Q 



Q 2 + Ml 



Q 2 /3 



(132) 

There are four distinct T 4 In T contributions to ^2k F ■ 
The first one comes from II 2 , and and is obtained by form- 
ing a cross-product of the two terms in Eq. l|131fl • In the 
limit of \Q\ > \Cl m \, 



n„ 



16 



[2gin|Q|+7r|fi m |] , 



(133) 



and this combination results in a nonanalytic term in 

^2k F 



^2k F OC 



TV|O m | 3 / ^ocT 4 ln^. 
^ J\n m \ Q T 



(134) 



Contrary to the q = contribution, however, the explicit 
T— dependence of II [entering via the second and third 
term in Ea. H13U| ] give three more T InT contributions 
in the thermodynamic potential. The second and third 
contributions come from cross-products of the T 2 and T 4 
terms in Ea. (|129[l with IIo. The logarithmic terms here 
come, respectively, from JdQ/Q — > ln.fl m — > InT and 
from T^2 n l/|f2 TO | — > InT. The fourth contribution 
comes from the square of the T 2 term in Eq. (|129f) . Each 
of these four contributions is of order T 4 In T, such that 



^2kf 



Using 



^{u F U(2k F )fNT(^- 
8 \vfKf 



3 E 4 



i=l 

(135) 



T 3 ^|a„|ln|a„|-^T 4 ln^, 



1 1 , 



T 4 J2^ L T 4 In 



Ep 
T ' 



we find 

7i : 
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14 2 
45' 74 ^9- 



45' ' 9 
Collecting all contributions, we finally obtain 



4 / rp 

S 2feF = — (v F U(2k F )) 2 NT — — 

5 \VFkF 



111 ■ 



Ep 
T 



(136) 



(137) 



(138) 



It is instructive to compare this result with small q contri- 
bution to S. We recall that II (e m , q) for \e m \ <C q <C kp 
is given by Ea. Hl(Jl|l . The logarithmic singularity in 5 is 
obtained by taking the cross-product of the |e m | /q and 
(|e m | /qf terms in Ea. fTT)T)l 



S cx TJ2 km I 3 [dqq 2 /q 3 OC T 4 In 
Calculating the prefactor, we obtain 



Ep 
T ' 



S ?=o = y (v F U{Q)) 2 NT 



^-Vln^ 

VpkF J T 



This almost coincides with Eq. i|138|) , the only difference 
is between the prefactors C/ 2 (0) and U 2 (2kp). 
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